Charge order at the frontier between the molecular and solid states in Ba3NaRu209 
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We show that the valence electrons of Ba3NaRu209, which has a quasi-molecular structure, 
completely crystallize below 210 K. Using an extended Hubbard model, we show that the charge 
ordering instability results from long-range Coulomb interactions. However, orbital ordering, metal- 
metal bonding and formation of a partial spin gap enforce the magnitude of the charge separation. 
The striped charge order and frustrated hep lattice of RU2O9 dimers lead to competition with 
a quasi-degenerate charge-melted phase under photo-excitation at low temperature. Our results 
establish a broad class of simple metal oxides as models for emergent phenomena at the border 
between the molecular and solid states. 

PACS numbers: 75.25.Dk,75.47.Lx 



Strong many body correlations between electrons in 
solids result in extremely diverse properties including 
insulating antiferromagnetism, charge order and super- 
conductivity. This is especially true in organic based 
materials like fullerides or the k-{BEDT-TTF)2X or 
(Pd[dmit]2)2 salts, due to the narrow electronic band- 
width resulting from the overlap of molecular orbitals [I]- 
[5]. However, the contacts between molecules are weak, 
and organic bonding allows low-energy rotational and 
conformational degrees of freedom. Poorly crystalline or 
complex, low symmetry structures are thus com- 

mon. This problem can be compounded by extreme air 
sensitivity, the invisibility of hydrogen to X-ray diffrac- 
tion, and the incoherent background which it generates 
in neutron scattering experiments. Consequently, nearly 
all theoretical models of these materials rely on approx- 
imations that treat each molecular unit as effective sites 
on high symmetry lattices [31 [7] . These models predict 
that Coulomb interactions drive electron localisation, 
and that geometrical frustration is responsible for the ob- 
served competition between ground states [S]. However, 
the material-specific interplay between electronic corre- 
lations, local-degrees of freedom, and the lattice remains 
a significant challenge [S]- 

fn this Letter, we show that simple metal oxides can 
provide new insight into the electronic properties which 
emerge between the molecular and solid states. Our 
model compound is Ba3NaRu2 09, which has the 6H 
hexagonal perovskite structure (Figs, la and lb) contain- 
ing isolated Ru2'^"'"09 dimers of face-sharing octahedra 
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FIG. 1: (color online) a) CO P2/c structure of Ba3NaRu209 
as determined from synchrotron powder X-ray diffraction pro- 
jected down [100], Sodium atoms are yellow spheres. Barium 
atoms are green. The symmetry inequivalent RU2O9 dimers 
are shown in red (Ru2'''09) and in blue (Ru2^09). (b) P2/c 
structure projected down [001] showing the frustrated trian- 
gular layers, (c) Results of magnetic susceptibility measure- 
ments for Ba3NaRu209 in a 1000 Oe field, showing fit to Ru''''' 
dimer correlations below 210 K. The inset (same units) shows 
the hysteresis on warming/cooling. 



which form triangular layers [ID] . Previous reports have 
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indicated poorly-metallic conduction around room tem- 
perature and a metal-insulator transition [llj associated 
with a structural distortion below 210 K. We synthesised 
single crystal and powder samples of Ba3NaRu20g which 
were characterised by magnetisation and high resolution 
synchrotron X-ray and time-of-flight neutron diffraction 
measurements. We confirmed that the room temperature 
structure has hexagonal P63 /mmc symmetry in which all 
Ru sites are equivalent. The magnetic susceptibility was 
observed to approximately halve on cooling through a 
first-order transition at Tco = 210 K, consistent with the 
opening of a charge or spin gap (Fig. fc). This is accom- 
panied by the formation of an apparently C-centred, or- 
thorhombic 2.a x 2.^3a x c superstructure |10j . however 
our high resolution X-ray synchrotron powder diffrac- 
tion measurements, performed using ID31 at the ESRF, 
detected a previously unreported monoclinic distortion 
(Fig. 2a) and further weak reflections that showed that 
the true space group is P2/c. The structural model was 
refined using X-ray data from a twinned single-crystal, 
and full details are given in the supplemental material. 
The key feature of our low temperature P2/c structure 
is a splitting of the Ruj "^^Og dimers into two symmetry 
inequivalent units, which suggests charge ordering (CO) 
(Fig. lb). Motivated by experiments on true molecular 
materials, we investigated the driving force for this tran- 
sition using the following extended Hubbard model, with 
the dimers as effective sites. 
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Here t is the intcrdimer hopping integral, C/g/ / is the 
on-site (dimer) Coulomb repulsion and V is the inter- 
dimer Coulomb repulsion. In contrast to Mott insulating 
organic materials such as the k-{BEDT-TTF)2X family, 
where V is commonly neglected [9], the experimentally 
determined stripe order could only be reproduced when 
V > Ueff. This implies that the formation of molecular 
orbitals results in a reduced Uef / [7] , and that the inter- 
site electron-electron interactions play a fundamental role 
in the ordering mechanism. However, without breaking 
the three-fold lattice symmetry, a striped groundstate 
is still highly degenerate. We discovered that the full 
three-dimensional P2/c structure could be rationalised 
by considering the tetrahedral linkages between the hep 
network of dimers. Each tetrahedron of dimers contains 
two Ru2'''09 and two Rug^Og units (SI). Hence, so- 
called 'ice rules', which minimise the electrostatic energy 
of the molecular building blocks [TH [T3] , act as directors 
for the fluctuating stripe order stabilised by the reduced 
Ueff- On the macroscopic scale, the twinning in single 
crystal samples has exactly the degeneracy predicted by 
this picture, as six orientational stripe domains of equal 
population are found. Our observations thus establish 
the Ru20g dimers as the physically meaningful units in 
Ba3NaRu20g. This analysis makes no prediction of the 
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FIG. 2; (color online) (a) Splitting and peak broadening ob- 
served in Ba3NaRu2 09 by synchrotron powder X-ray diffrac- 
tion upon cooling into the P2/c CO phase, (b) The Ru-Ru 
distance has a linear dependence on oxidation state in the 
Ba3ARu209 family of compounds with A = Na"*", Ca^^, 
Nd^^ and Ce''"'", which have similar ionic radii. The for- 
mal ruthenium oxidation states are 5.5^, S"*", 4.5^ and 4+, 
open symbols are from the P2/c phase determined in this 
work. The shaded region shows the standard deviation of the 
calibration fit. (c) A similar linear dependence upon formal 
oxidation state is found for the average Ru-O bond distance 
in the Ba3ARu209 materials. The symbols are as above, 
(d) Coordination environment in the two inequivalent RU2O9 
dimers. 

orbital ordering. 



The Ru®^ octahedra are compressed implying t^g 



size of the charge separation between the sites. How- 
ever, our precise refinement of P2/c structure allowed us 
to examine structural parameters on the single-ion and 
molecular levels, and compare them to other members of 
the Ru20g family of known oxidation states. The valence 
of each individual Ru cation is reflected [1^ by a linearly 
varying <Ru-0> distance (Fig. 2b), and the refined val- 
ues for the P2/c structure are in agreement with 'inte- 
ger' separation into 5"'" and 6^ oxidation states. Further- 
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more, detailed examination of the coordination environ- 
ment (Fig. 2d) of the two sites provides evidence for or- 
dering of orbital degrees of freedom. In order to separate 
out the effect of the intrinsic trigonal site distortion from 
electronic effects, we calculated the average length of per- 
pendicular bond pairs. For the Ru^"*" sites, which are or- 
bitally non-degenerate with a configuration, these are 
extremely regular, and lie in the range 1.944(5)-1. 9815(8) 
A. In contrast, the t^g Ru^"*" sites have four long bonds 
with average lengths in the range 1.943(5)-1. 991(5) A 
and two short bonds of average length 1.866(5) A. An 
orbital ordering distortion therefore lifts the degeneracy 
of the doubly occupied t2g orbitals. In addition to these 
single-ion indicators of integer charge order, the inter- 
dimer Ru-Ru distances also vary linearly with oxidation 
state (Fig. 2c). This implies a well defined 'bond order' 
determined by the occupation of molecular orbital states 
and allows a precise estimate of the charge separation 
into 4.86(16)+ and 6.11(16)-h states. The Ru-Ru bond 
in the Ru2^0g dimers is as short as that found in Ru 
metal (2.62 A). Metallic bonding would make these sites 
non-magnetic, which is in accordance with the magnetic 
susceptibility below 210 K, which could be modelled as- 
suming only S=l correlations from the Ru2^0g dimers 
(Fig. Ic). Finally, we used the precise P2/c coordinates 
to perform LDA-I-U calculations with U = 3 eV, which 
confirm that complete disproportionation into Ru2'''09 
and Ru2'''09 pairs is found (Fig. 3a). This is stabilised 
by the S = dimer states, as antiferromagnetic coupling 
gives by far the lowest energy groundstate. A small gap 
at the Fermi energy is found in accordance with Ref . [TT] . 
These calculations confirm that the simplest molecular 
constraint, the desire to form closed shell configurations, 
stabilises the complete charge segregation found below 
210 K . We emphasise that our experimental and theo- 
retical quantification of the degree of charge localisation 
in Ba3NaRu2 09 would have been effectively impossible 
for organic materials. 

The integer charge separation motivates comparison with 
frustrated Ising spin models on the hep lattice, which 
have the same groundstate configuration [Tl] . Geometri- 
cal frustration of electronic degrees of freedom in molec- 
ular solids is additionally predicted to result in com- 
petition between charge ordered and disordered states 
[Ulin]. Although the P2/c phase of Ba3NaRu2 09 has 
large (-0.1 A) atomic displacements, we discovered that 
the charge order is unusually sensitive to external per- 
turbations. A comparison between the lattice parame- 
ters extracted from X-ray (ID31 at ESRF) and neutron 
(WISH at ISIS) diffraction experiments is shown in Fig. 
3b. Below —40 K these deviate strongly from each other. 
This was investigated by cooling a sample to 10 K before 
continuously measuring diffraction profiles with 31 keV 
X-rays. We found that the P2/c charge ordered struc- 
ture transforms continuously to a higher symmetry C2 / c 
structure in which only one Ru site is present. This phase 
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FIG. 3: (color online) (a) Site specific ruthenium density of 
states determined by LDA+U calculations for the charge or- 
dered P2/c phase of Ba3NaRu2 09, a small insulating gap is 
stabilised, in accordance with Ref. (12). (b) Temperature 
dependence of lattice parameters determined by Rietveld re- 
finement of powder synchrotron X-ray and neutron powder 
diffraction. Lines show fits used to extrapolate T = data 
points used in Fig. 4 



was found to be stable for periods of several hours at 10 
K, during and after irradiation, and was not observed 
by neutron powder diffraction down to 1.6 K, which also 
showed no evidence for long-range magnetic order. Band 
structure calculations predict metallic conduction in this 
phase. 

The synchrotron X-ray data collected during the melt- 
ing process show that local and lattice degrees of free- 
dom are inextricably linked. As shown in Fig. 4a, the 
average Ru^+/Ru^+ charge order, quantified by the dif- 
ference between the Ru-Ru distances in the two dimers, 
decreases rapidly with an exponential time constant of 
75(8) s. After 120-180 s, this charge order and the cor- 
responding {h + k) = odd superstructure reflections are 
completely suppressed showing that the average struc- 
ture has transformed to C2/c symmetry. The very high 



4 



resolution of the synchrotron diffraction data makes a 
transition to an intermediate charge glass of localised but 
spatially disordered charges unlikely [16], as the micro- 
scopic lattice strain, which measures local variations in 
lattice parameters, also rapidly vanishes. The lattice dis- 
tortion which breaks frustration, parameterised by the 
h/ a ratio, is coupled to the local order, and relaxes on a 
similar timescale (88(3) s). 

A possible explanation for the CO melting is the dimen- 
sional reduction into stripes |17| patterned by the ice- 
rule constraints. X-rays impinging on the sample gen- 
erate a low concentration of holes by the photoelectric 
effect. Based on our minimal electrostatic model, these 
may freely delocalise along the stripe direction, provided 
that t > y — melting the charge order. That this 

charge melted phase exists in delicate balance with the 
charge order, is shown by measurements on warming, as 
the P2 1 c phase is recovered on warming above 40 K with 
full recovery of the lattice micro and macrostrains. Fu- 
ture work will be needed to determine the factors which 
set this temperature scale. However, we note that x{T) 
also reproducibly shows a glassy transition |10j at similar 
temperatures (Fig. Ic), which might indicate an intrinsic 
instability of the CO, even in the absence of irradiation. 
Unlike other charge-ordered metal oxides, Ba3NaRu2 0g 
shows both fully crystallized and liquid-like valence elec- 
tron states at low temperatures. However, the classic 
examples of oxide charge order are all found in struc- 
ture types which have infinite lattices of corner, edge or 
face sharing coordination polyhedra, like perovskites or 
spinels. The mechanism for charge order is somewhat 
controversial in these materials, as only fractional charge 
separation is found [TH], with strong evidence for Fermi 
surface nesting [12]. In contrast, our results shows that 
long-range electrostatic interactions are responsible for 
charge ordering in Ba3NaRu209. Additionally, the com- 
paratively simple crystal structures highlight the role 
played by single-ion and molecular degrees of freedom. 
Our results thus lend weight to a recent theoretical in- 
vestigation of K-ET2Cu2(CN)3, which showed that intra- 
molecular degrees of freedom must be important in this 
related class of materials [2D]. 

In summary, our results show that a combination of 
single-ion and molecular degrees of freedom help to sta- 
bilise charge order in Ba3NaRu2 09. Our results are remi- 
niscent of the true (organic) molecular solids and provide 
a model system for studying frustrated strongly corre- 
lated electrons. We note the presence of a vast library 
of unexplored related compounds [21', and predict the 
future discovery of strongly fluctuating spin-liquid or su- 
perconducting groundstates. 
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ments and the HZB, EPSRC and Leverhulme trust for 
support. We thank S. Parsons for assistance with the 
single crystal data collection, and CD. Ling, R. Valenti, 
V. Honkimaki and T. Chatterji for useful discussions. 
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FIG. 4: (color online) (a) Difference between the refined Ru- 
Ru distance for the two dimer types in Ba3NaRu209as a func- 
tion of irradiation time at 10 K. The {h + k)=odd superstruc- 
ture reflections (inset) also disappear, showing that the CO 
melts. Results of fitting exponential relaxation functions are 
shown; (b) Time dependence of the lattice macrostrain (b/a) 
ratio. The inset shows the immediate decrease in the lattice 
microstrain, determined from the tanS dependent contribu- 
tion to the Lorentzian peak shape. 
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FIG. 5: (color online) SI The hep lattice of RU2O9 dimers in Ba3NaRu209 has tetrahedral linkages between dimers in the 
stacking (c) direction. In the charge ordered phase, two of each charge state are found in each tetrahedral unit, which minimises 
the electrostatic energy of the charge ordered state. 
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TABLE I: Details of single crystal X-ray data collection and refinement of Ba3NaRu2 09 at 110 K 



Formula 
Cell setting, Space group 
Temperature(K) 
a,b,c (A) 

/3 

Volume (A^) 
Dc (mg.m-^) 
Radiation type 
No. reflections for cell parameters 
H (mm~^) 
Crystal form, colour 
Crystal size (mm) 
Diffractometer 
Data collection method 
Absorption correction 



No. of measured, independent and observed refiections 

ft (°) 

Rint 

Range of h,k,l 

Refined on 

R, wR(F), S 

Cutoff: I > 
No. of reflections 
No. of parameters 

Weighting scheme 



14.48497(6) 



Ba3NaRu209 
Monoclinic, P12/cl 
110 K (2) 
5.84001(2) , 10.22197(4) 
90.2627(3) 
864.69(2) 

6.00 
Mo-Ka 

725 
16.948 
hexagonal plate, black 
0.13 X 0.45 X 0.45 
Bruker Apex2 

w + <p 
Spherical 
0.0616 
0.1174 
11855, 2549 
30.5293 
0.067 
-8 h 8\ 
fc 14 
/ 20/ 
F 

0.028, 0.032, 1.01 
3.00o-(I) 
2651 
80 

1 - (AFob./6 X AF,,tf] 



21 2 



(A/a)max 

Ap 

max ; Apmin (e A ' ) 



w — W X 

w' = [PoTiix) + PiTiix) + . . . P„-iT;,_,{x)]- 

where Pi are the coefficients of a 
Chebychev series in ti{x), and x = Fcaic/-Pcaicn 
0.0010 
1.97, -1.36 
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TABLE II: Twin laws and refined twin fractions 



Twin law Twin fraction 




0.1656(15) 
0.165(3) 
0.168(3) 

0.1674(15) 
0.169(3) 
0.169(3) 



TABLE III: Refined atomic coordinates and displacement parameters for Ba3NaRu209 at 110 K from single crystal X-ray 
diffraction. Anisotropic displacement parameters were only refined for the barium atoms. 



Atom 


X 


y 


z 


^iso/ equiv 


Na(l) 











0.0044(5) 


Na(2) 


0.5 


0.5 





0.0052(5) 


Ba(l) 





0.00689(4) 


0.25 


0.0051 


Ba(2) 


0.5 


0.48195(5) 


0.25 


0.0052 


Ba(3) 


0.50537(14) 


0.15771(3) 


0.91817(2) 


0.0053 


Ba(4) 


0.00413(14) 


0.68674(4) 


0.90174(2) 


0.0057 


Ru(l) 


0.49931(19) 


0.16087(3) 


0.15941(2) 


0.00395(9) 


Ru(2) 


0.00241(9) 


0.65572(4) 


0.15151(3) 


0.00426(9) 


0(1) 


0.7317(6) 


0.2404(4) 


0.2461(6) 


0.0054(8) 


0(2) 


0.2065(7) 


0.7367(4) 


0.2361(3) 


0.0040(8) 


0(3) 





0.4843(6) 


0.75 


0.0095(12) 


0(4) 


0.5 


-0.0109(6) 


0.75 


0.0096(12) 


0(5) 


0.2694(14) 


0.0793(9) 


0.4087(5) 


0.0128(19) 


0(6) 


0.8090(9) 


0.5704(5) 


0.4236(3) 


0.0083(9) 


0(7) 


0.0074(15) 


0.1905(4) 


0.9177(3) 


0.0045(7) 


0(8) 


0.4920(16) 


0.6846(4) 


0.9122(3) 


0.0040(7) 


0(9) 


0.7365(10) 


0.0857(7) 


0.4108(4) 


0.0025(14) 


0(10) 


0.2827(8) 


0.6028(5) 


0.3892(3) 


0.0080(9) 



